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ABSTRACT 

Carbonaceous chondrite meteorites are known for having high water and organic material contents, 
including amino acids. Here we address the origin of amino acids in the warm interiors of their 
parent bodies (planetesimals) within a few million years of their formation, and connect this with the 
astrochemistry of their natal protostellar disks. We compute both the total amino acid abundance 
pattern as well as the relative frequencies of amino acids within the CM2 (e.g. Murchison) and 
CR2 chondrite subclasses based on Strecker reactions within these bodies. We match the relative 
frequencies to well within an order of magnitude among both CM2 and CR2 meteorites for parent body 
temperatures < 200°C. These temperatures agree with 3D models of young planetesimal interiors. 

We find theoretical abundances of approximately 7x10® parts-per-billion (ppb), which is in agreement 
with the average observed abundance in CR2 meteorites of 4±7xl0®, but an order of magnitude higher 
than the average observed abundance in CM2 meteorites of 2 ± 2x10"^. We find that the production of 
hydroxy acids could be favoured over the production of amino acids within certain meteorite parent 
bodies (e.g. CII, CM2) but not others (e.g. CR2). This could be due to the relatively lower NH 3 
abundances within CII and CM2 meteorite parent bodies, which leads to less amino acid synthesis. 

We also find that the water content in planetesimals is likely to be the main cause of variance between 
carbonaceous chondrites of the same subclass. We propose that amino acid abundances are primarily 
dependent on the ammonia and water content of planetesimals that are formed in chemically distinct 
regions within their natal protostellar disks. 

Key Words: astrobiology — astrochemistry — meteorites, meteors, meteoroids — molecular data 
— methods: analytical 


1. INTRODUCTION 

One of the central questions about the origin of life 
is whether the delivery of exogenous water and organ¬ 
ics to planets may have ultimately enabled it. Car¬ 
bonaceous meteorites certainly deliver a rich mixture of 
biomolecules to the Earth today. These organic mate¬ 
rials, including amino acids, could have formed in a va¬ 
riety of environments such as the interstellar medium, 
protostellar disks out of which planets and stars form, 
and within the planetesimals and asteroids that built the 
rocky planets. It is well known that these extraterrestrial 
sources could have contributed a very significant fraction 
of the organics present on the young E arth’s surface at 
the dawn of life ( Chyba fc Sagan|[l992 ). 

A direct connection that we have with the conditions 
that created such early amino acids comes from their 
measured abundances and relative abundance patterns 
within meteorites. A good example is the extensive 
chemical analysis that has been done of carbonaceous 
chondrite meteorites such as the famous Murchison me¬ 


teorite (Cronin & Pizzarello 19831. The parent bodies 
of meteorites namely planetesimals which are typically 
up to 50 - 100 km in size, were formed during the cre¬ 
ation of the solar system. Their composition consists 
of rock and ice together with various simple molecules 


^ Origins Institute, McMaster University, ABB 241, 1280 Main 
St, Hamilton, ON, L8S 4M1, Canada 

^ Department of Physics and Astronomy, McMaster Univer¬ 
sity, ABB 241, 1280 Main St, Hamilton, ON, L8S 4M1, Canada 
^ E-mail: alyssacobbl07@gmail.com 
^ E-mail: pudritz@physics.mcmaster.ca 
® E-mail: ben.pearce@alumni.ubc.ca 


(eg. HCN and ammonia) and organics (eg. formalde¬ 
hyde) that were inherited from the icy, dusty materials 
in disks and interstellar gas. Such bodies also contained 
various radionuclides whose decay heated their interiors 
up to 150°C for several million years. Given time and 
meteorite impacts, their organic content has been con¬ 


tinually introduced to the planet (Chyba fc SaKan||I992 


Sephton 2002 Botta & Bada 2002'Clavin et al. 2011). 


i'his abundance pattern may have provided important 
constraints on the nature of the genetic code that ulti- 


mately appeared ( 

Higgs & Pudritz 

2009 

I- 

I’he lormation ot complex molecu 

les IS still not well un 


derstood. Experimental work by Lerner et al. (1993) has 
shown that carbonaceous chondrites like Murchison (a 
so-called CM2-type) may be explained by a Strecker-type 
synthesis wherein aldehyde molecules (such as formalde¬ 
hyde) in the presence of water, ammonia and HCN 


give rise to amino acids (such as glycine). Lerner et 
al. (|1993 ) shows that Strecker synthesis retains the ob¬ 


served highly deuterated amino acids found in Murchi¬ 
son, whose source ultimately is the D-enriched precur¬ 
sors in the interstellar medium. Moreover, with few ex¬ 
ceptions, it is found that the major amino acid that is 
produced in any reaction is the one expected from the 
associated aldehyde. Not all classes of amino acids are 
necessarily formed in such reactions however, and there 
is some evidence that the so-called n-w species may have 
been formed by catalytic processes associate d with so- 


cal led Fi scher-Tropsch-Type (FTT) reactions (Burton et 


al. 


2012 ). 


he observation of planetesimals and complex organic 
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molecules in protostellar disks is now on the cusp of a 
revolution in our understanding with the advent of the 
ALMA observatory. As one example, it will enable obser¬ 
vations of planetesimal systems formed in disks around 
other stars as well as the capability to detect carbon¬ 


beari ng species in the later debris disks (Dent et al. 


20141. ALMA has also been able to detect glycolaldehyde 


(the precursor to the amino acid serine) in abundances 
10-15 times less than methyl formate and 2-3 times more 
than ethylene glycol in the warm gas surrounding the 


proto stellar binary IRAS 16293-2422 (Jprgensen et al. 


20121. Such aldehydes would be incorporated into plan- 


etesimals and comets where they could play a crucial role 
in amino acid synthesis, as we will show. 

The overarching goal of this work is to thermodynami¬ 
cally model the synthesis of amino acids inside meteoritic 
parent bodies and to begi n to link this to the astr ochem- 
istry of their natal disks. Cobb & Pudritz (2014) (here¬ 
after Paper I) took a closer look at overall amino acid 
abundances and the relative frequencies of amino acids 
within specific meteorite subclasses. This comprehensive 
view of amino acid abundance data provided us with the 
constraints that we now use in this computational study 
to reproduce these observed abundance patterns. 

Early theoretical work on the synthesis of a mino acids 
in hydrother mal environments was done by [Schulte 


Shock ([1995 ). Using principles of thermodynamics and 


Cibbs free energies, they showed mathematically that 
the synthesis of amino acids via a Strecker synthesis-like 
pathway is energetically fa v orable in hydrothermal con¬ 
ditions. Amend & Shock (19981 used thermodynamic 
considerations and Gibbs free energies to produce a the¬ 
oretical yield of amino acids in hydrothermal ecosystems. 
They found that the synthesis of 20 proteinogenic amino 
acids is energetically favorable in hot (100°C) hydrother¬ 
mal environments, as op posed to their synthesis at sur¬ 
face seawater conditions. Schulte & Shockj (|2004|) shifted 
focus from hydrothermal systems to meteoritic parent 
bodies. Using theoretical compositions of carbonaceous 
chondrite parent bodies, they calculated that various or¬ 
ganic materials present in carbonaceous chondrites could 
have been synthesized during aqueous alteration of the 
parent body. 

In this work, we investigate the theoretical synthesis 
of proteinogenic a-amino acids via Strecker-type path¬ 
ways inside meteoritic parent bodies. These bodies are 
important for amino synthesis since their interiors are 
heated to temperatures exceeding 100° C for more than 
1 Myr due to radioactive heating, as shown by full 3D 


hydrodynamics simulations of their interiors (Travis & 


Schubert 20051. Moreover, the time scale for these reac¬ 


tions which synthes ize amino acids is no longer than a 
few thousand years (Peltzer et al.|1984). Therefore it is a 
good assumption that these processes will have achieved 
thermodynamic equilibrium, making the Gibbs free en¬ 
ergy minimization a desirable approach. Our model in¬ 
herently assumes this chemical equilibrium, and outputs 
a series of amino acid concentrations which minimize the 
Gibbs free energy. We model parent body interiors using 
tempe ratures and pressures given in [Travis fc Schubertj 
(2005) as well as cometary abundance data of organics 
for the initial concentrations in our model. 

We compare our modeled amino acid abundances 
to observed total abundances in different subclasses of 


carbonaceous chondrite, measured in parts-per-billion 
(ppb). Additionally, in Paper I we calculated the av¬ 
erage relative frequencies of various amino acids normal¬ 
ized to glycine, measured in moles of amino acid/moles 
of glycine. The amino acid abundances measured in hy¬ 
drolyzed (instead of unhydrolyzed) hot water extracts 
were used for these calculations due to a larger quan¬ 
tity of available data. These average frequencies have 
since been refined and standard deviations have been 
added. We note that the variation in the measured rela¬ 
tive amino acid frequencies is greater than any discrep¬ 
ancies introduced from using amino acid abundances in 
hydrolyzed meteorite extracts. We calculate the pattern 
of relative frequencies in theoretical yield, and match it 
to observation in two carbonaceous subclasses, CM2 and 
CR2. 

In this work, we limit our investigation to proteino¬ 
genic a-amino acid abundances. This is primarily be¬ 
cause the thermochemical database with which we work 
currently only contains data for the 20 proteinogenic 
amino acids encoded by the universal genetic code— 
with the exception of one non-proteinogenic a-amino 
acid. It is important to note that there are apprecia¬ 
ble quantities of non-proteinogenic amino acids (such 
as amino-isobutyric acid, isovaline and a-aminobutyric 
acid) and non-a-amino acids (such as /3-alanine) detected 
in carbonaceous chondrites. In Appendix I, we analyze 
one case: a-aminobutyric acid. In addition, we cur¬ 
rently only consider amino acids which are synthesized 
via Strecker-type reactions involving aldehydes. Amino 
acids such as amino-isobutyric acid and isovaline, while 
both highly abundant in carbonaceous chondrites, are 
synthesized via Strecker-type reactions with ketone pre¬ 
cursors. 

Our results give very good agreement with the ob¬ 
served data base of amino acid frequencies. In trying 
to understand the absolute abundances of amino acids 
in the various meteorite classes predicted in our models, 
we propose that the variations in the water and ammonia 
content of planetesimals formed at varying radii within 
natal protoplanetary disks would naturally account for 
these differences. This is a new picture that is quite dif¬ 
ferent than the current ’’onion” structure of differentiated 
chemistry within a planetesimal. 

The layout of this paper is to first quickly review our 
observational results from Paper I, summarize the the¬ 
oretical and computational methods that we used, and 
then to discuss our results and comparison to the data 
for both the abundances and relative frequencies of our 
models. Finally, we discuss the results in a larger frame¬ 
work and introduce our new picture in the last section of 
the paper. 

2. CO]V[PILATION OF METEORITIC AMINO ACID DATA 

Paper I presented an extensive collation of available 
meteoritic data on amino acid abundances and relative 
frequencies in carbonaceous chondrites. Carbonaceous 
chondrites are the class of meteorite that contain the 
highest content of both water and organic materials, in¬ 
cluding amino acids. These abundances are seen to vary 
systematically with meteoritic type. In this section we 
give a quick summary of meteorite classification, followed 
by a summary of their amino acid abundances (the reader 
is encouraged to read Paper I for the details). 



























Simulating the Synthesis of Amino Acids in Meteorite Parent Bodies 


3 


2 .1. Meteorite Classification - a Quick Guide 

Meteorites are assigned a number 1-6 indicating petro¬ 
logic type, in addition to a subclass designation Cl, CM, 
CR, CV, CO, CH, CB, and CK. Petrologic type carries 
with it information regarding degrees of internal chem¬ 
ical alteration undergone by a meteoritic parent body. 
Petrologic type 1 chondrites are the most aqueously al¬ 
tered. This means the liquid water within the parent 
bodies has chemically altered their makeup. Petrologic 
type 2 meteorites have undergone significant degrees of 
aqueous alteration, though not to the extent of type 1 
meteorites. Petrologic type 3 meteorites most closely re¬ 
semble the solar nebulae in which they formed; they have 
undergone relatively little to no degrees of alteration. All 
carbonaceous chondrites subclasses from petrologic types 
1-3 are known to contain high concentrations of water 
and organics. Petrologic types 4-6 indicate increasing 
degrees of thermal metamorphism. We restrict our mod¬ 
eling to meteorites of type 1-3, as samples in this regime 
have been aqueously altered and produced amino acids 
via Strecker synthesis. 

Cl meteorites, all petrologic type 1, are among the 
most aqueously altered carbonaceous chondrites. They 
are generally associated with a high er range in temper ¬ 
atures, 50°C to 150°C , according tolKrot et al.| (|2006|), 
Zolensky^^etaL (|1993 ), Huss et al. (2006), and Grirnm 


&: Mcbween 1 TD8P ), and they rank approximately fourth 


highest tor average total amino acid abundances by car¬ 
bonaceous chondrite subclass. 

The CM and CR subclasses express the highest to¬ 
tal amino acid concentrations, greater by several orders 
of magnitude than any other subclass. They also ex¬ 
press the greatest variety of proteinogenic amino acids. 
These subclasses generally correspond to petrologic type 
2 , with several exceptions, and are associated with alter¬ 
ation temperatures in the broad range of 0°C to 240° C. 


some studies dZolensky et al.|1993| 

Huss et al. 12006 

Krot 

et al. 

2UU6 Grimm & iVlcbween[ 

1989 

), but one study 


of 240° C (Busemann et al. 2007). Peak metamorphic 
temperatures of chondrites are typically estimated from 
phase c hanges or reactio n s that occur at specific t emper- 
atures (Huss et al. ]2006l ). [Busemann et al.|(|2007[) on the 


other hand use the correlation between the D band width 
(obtained from Raman spectroscopy) and peak metamor¬ 
phic temperature to base their estimates. CR meteorites, 
spanning petrologic types 1 -3, also have a discrepan cy in 
classification temperatures. Busemann et al. ( 2007|) give 


a peak metamorphic temperature of 240°C, [Huss et al.| 


(200_6) cite an upper temperature bound of 150-300°C, 


and Zolensky et ah ( 1 993^ s tate a range of 50°C to 150°C. 


iVlorlok & Libourei|| (|2U13|) provide a broader range in 
classihcation temperatures, 25°C and 300°C, associated 
with “all stages of aqueous alteration.” 

CV, CO, CH and CB meteorites are all petrologic type 
3. These meteorites are associated with peak metamor¬ 


phic temperatures ranging fro m 200-60 0°C (Busemann 


et al. 


T 


2007 


I 


Huss et al. 2006), though Zolensky et al. 


1993) has CV3 peak teiriperatures at a niuch lower 50°C. 


hese chondrites are thought to be the most pristine 
samples of material we have from the time of the solar 
system formation. They are relatively unaltered, either 


via aqueous alteration or thermal metamorphism—which 


occurs at hif 

her temperatures of perhaps ~ 500° — 950° C 

(Huss et al. 

2006 

) with no water present. With the ex- 

ception of C 

H meteorites which are the third most abun- 


dant subclass in total amino acids, we observe low levels 
of amino acids in these meteorites. 

Lastly, the CK meteorites span petrologic types 3-6 
and are associated with temperatures from 250-600° C 
( Burton et al.|2015 ). A recent study has found one CK3 
meteorite to be c omparable in amino acid abundances 
to the CIl types (Burton et al. 2015). The CK4-6 me¬ 
teorites in the sarne study were found to have very low 
abundances of amino acids. 


2 .2. Amino Acid Abundances 


In Figure 8 of Paper 1, we plot average amino acid 
abundances in parts-per-billion (ppb) across different 
carbonaceous chondrite subclasses. The histogram style 
abundance plot shows average amino acid abundance for 
various proteinogenic amino acids for each of the six sub¬ 
classes of carbonaceous chondrite we consider. This al¬ 
lows us to identify graphically in which subclasses there 
is a greater concentration of amino acids. The C03 and 
CMl/CRl sections have average abundances on the or¬ 
der of 10^ ppb. The average abundance of each amino 
acid for CV3 and CIl subclasses is on the order of 10^ 
ppb. CM2-types contain 10^ ppb per amino acid on av¬ 
erage. The CR2-types show the greatest variation in 
abundance, each average ranging from 10^-10® ppb. The 
CH3, CB3 and CK3 types were not included in Figure 
8 of paper I, but have since been found to contain 10^- 
10 ^ ppb, 10^-10^ pp b and 10^ ppb per amino acid on 
average, respectively ([Burton et al.||2013[ |2015|). From 
these two studies and I'igure 8 of Paper 1, we identified 
the petrologic type 2 meteorites, both CM2 and CR2, as 
containing greater total abundances of amino acids than 
either the petrologic type 1 or type 3 meteorites. This 
is in agre ement with observations and conclusions b y, for 
example, Martins et al. (2007b), Glavin et al. (2011), and 
Burton et al. 12012). 


2.3. Amino Acid Frequencies 


In addition to modeling the total amino acid yield, we 
investigate the relative frequencies of amino acids. Fig¬ 
ures 9 and 10 in Paper I show the relative frequencies of 
13 amino acids relative to glycine in the CM2 and CR2 
subclasses, respectively. As the CM2 and CR2 subclasses 
contain the greatest concentrations and variety of amino 
acids, we selected these groups as a starting place for 
thermodynamic modeling. We normalize to glycine as 
this amino acid is classically one of the most abundant 
found in meteorites, it is the easiest to synthesize and is 
most energetically favorable, and it is achiral. In both 
subclasses, alanine is the second most abundant amino 
acid and glutamic acid is the third. Among the CM2 
chondrites, there are also significant quantities of serine, 
proline, lysine, threonine, valine, isoleucine, leucine, as¬ 
partic acid, tyrosine, and phenylalanine. The CR2 sam¬ 
ples show great variety in frequency; valine, isoleucine, 
leucine, proline, and phenylalanine are fairly prevalent, 
and there are small quantities of aspartic acid, threonine, 
serine, and tyrosine. We use the same theoretical model 
used for calculating total amino acid yield to measure 
relative frequencies. 
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Figure 1. General form of a Strecker-type synthesis reaction. An aldehyde (R-CHO) reacts with hydrogen cyanide (HCN) and ammonia 
(NH 3 ) in aqueous solution (H 2 O) to produce an amino acid (H 2 N-RCH-COOH) and a new ammonia molecule. R denotes the side chain, 
which varies for different aldehydes/amino acids. 


3. THEORETICAL BACKGROUND 
3.1. Streaker Chemistry 

The general form of a Strecker-type reaction is shown 
by the mass-balance equation in Figure Following 
the experimental w ork discussed in the Introduction 
(Lerner et al. 1993 ), we simplify our equilibrium chem- 
istry scheme to focus on the relation between an alde¬ 
hyde and its main amino acid. A Strecker-type reaction 
is defined as the reaction between an aldehyde (R-CHO), 
hydrogen cyanide (HCN) and ammonia (NH3) in aque¬ 
ous solution (H2O) which produces an amino acid (H2N- 
RCH-COOH) and an ammonia molecule. There is an in¬ 
termediate step not depicted in Figure]^ in the Strecker 
synthesis of glycine, for example, reactants initially com¬ 
bine to create an aminoacetonitrile (NH2CH2CN), which 
then undergoes a one-way process called hydrolysis, re¬ 
sulting in the creation of glycine (C2H5N02). This gen¬ 
eral form may be used to define the underlying Strecker 
pathway for each of the proteinogenic amino acids laid 
out in Paper 1 . In each case, the side chain (R group) of 
a reactant aldehyde will be identical to the R group in 
the corresponding amino acid. 

This precise treatment of Strecker-type synthesis ne¬ 
cessitates the R-group equivalence of the reactant and 
the product. For each Strecker reaction producing a pro¬ 
teinogenic amino acid, we can balance the incoming atom 
count with the outgoing atom count, which allows for the 
definitive determination of the aldehyde required to cre¬ 
ate an amino acid. Table [T] shows this information for 
17 of the 20 proteinogenic amino acids: glycine (gly), 
alanine (ala), aspartic acid (asp), glutamic acid (glu), 
isoleucine (ile), leucine (leu), lysine (lys), serine (ser), 
threonine (thr), valine (val), tyrosine (tyr), phenylala¬ 
nine (phe), histidine (his), arginine (arg), tryptophan 
(trp), glutamine (gin), and asparagine (asn). We ex¬ 
clude proline (pro) due to its inconformity with the ba¬ 
sic Strecker-type treatment, as well as cysteine (cys) and 
methionine (met), which each contain a sulfur atom. For 
each listed amino acid. Table provides the amino acid 
name and molecular formula as well as the name and 
molecular formula of the requisite Strecker aldehyde. 

The information contained within Table may be used 
to construct the unique Strecker-type reactions resulting 
in the synthesis of each of 17 proteinogenic amino acids. 
Aldehyde and amino acid pairings may be input into the 
generic Strecker reaction shown in Figure [2 For example, 
the simplest Strecker pathway is for glycine, shown by 
Eq. below. 


CH2O + HCN 


NH, 


H2O 


H2N - HCH - COOH 


NH, 


( 1 ) 


Table 1 

Proteinogenic amino acids and their requisite Strecker 
aldehyde pairings. Includes the name and molecular 
formula of each amino acid and aldehyde. 


Amino Acid 

Aldehyde 

glycine 

formaldehyde 

C2H5NO2 

CH20 

alanine 

acetaldehyde 

C3H7NO2 

C2H40 

valine 

isobutyr aldehyde 

C5H11NO2 

C4H80 

leucine 

3 -methyl butanal 

C6H13NO2 

C5H10O 

isoleucine 

2 -methyl butanal 

C6H13NO2 

C5H10O 

serine 

glycolaldehyde 

C3H7NO3 

C2H4O2 

lysine 

5 -aminopentanal 

C6H14N2O2 

C5H11N0 

threonine 

2 -hydroxy-propanal 

C4H9NO3 

C3H6O2 

tyrosine 

p-hydroxyphenyl- 

C9H11NO3 

acetaldehyde C8H8O2 

phenylalanine 

phenyl-acetaldehyde 

C9H11NO2 

CgHgO 

aspartic acid 

3 -oxopropanoic acid 

C4H7NO4 

C3H4O3 

glutamic acid 

4 -oxobutanoic acid 

C5H9NO4 

C4H6O3 

histidine 

imidazole- 4 -acetaldehyde 

C6H9N3O2 

C5H6N20 

arginine 

N-( 4 -oxobutyl)guanidine 

C6H14N4O2 

C5H11N30 

tryptophan 

indol- 3 -ylacetaldehyde 

C11H12N2O2 

C10H9N0 

asparagine 

3 -oxopropanamide 

C 4 HgN 203 

C3H5NO2 

glutamine 

4 -oxobutanimide 

C5H10N2O3 

C4H7NO2 


Here the reactant is formaldehyde (CH2O) and the prod¬ 
uct amino acid is glycine (H2N-HCH-COOH, also seen in 
its collapsed form C2H5NO2 in Tablef^. The R group in 
both the aldehyde and the amino aciclis a single H atom. 
Note the traditional molecular structures shown in Fig¬ 
ure have been collapsed to their molecular formulae for 
ease of reading. 


3.2. Gibbs Free Energies 

Thermodynamically speaking, a chemical reaction is 
favorable if the Gibbs free energy of reaction is negative. 
Each Strecker-type reaction we consider has some Gibbs 
free energy of reaction associated with it (AGr)- In addi¬ 
tion to species-specific Gibbs energies of reaction, there is 
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a Gibbs free energy of formation (AG/) associated with 
each species involved in a reaction. AG,, may be com¬ 
puted at various temperature and pressure combinations 
using the AG/ for each species. These values are related 
via the relationship 


AGr = EG 


products 

/ 


_ ^^reactants 


( 2 ) 


The values of AG/ change with temperature and pres¬ 
sure for each species according to Eq. 

Gf(T,P) = a+bT+cTln{T)+dT^+eT^+f/T+gP (3) 


which provides us with the Gibbs coefficients a through g 
for each chemical species. Each term in this expression is 
derivable from the general definition of Gibbs free energy. 
Here, temperature (T) must be in Kelvin and pressure 
(P) in bars, and Gibbs energies return in Joules/mole. 
Eq. is the requisite Gibbs energy function a s req uired 
by our modeling software, discussed in Section [T^ 

4 . COMPUTATIONAL METHODS 



Figure 2. A plot of the reactant and product Gibbs energies of 
formation for the Strecker synthesis of glycine at 100 bar. The 
vertical grey dotted line represents the boiling point of water at 
100 bar. 


Our general strategy follows the outline described by 
Gibbs energies above. Initially we compute a series of 
Gibbs coefficients a-g for each chemical species involved 
in the Strecker-type synthesis of amino acids. Following 
this, we perform a minimization of the Gibbs free en¬ 
ergy for each individual reaction in the entire network 
of Strec ker e qua tion s. These two tasks are discussed in 
Sections |4.1| and |4.2[ respectively. 

4.1. Computing Gibbs Formation Energies 

We use a chemical thermodynamics database called 
CHNOSZ (version 0.9-9 (2013-01-01), authored by Jef¬ 
frey M. Dick, http://www.chnosz.net/). CHNOSZ con¬ 
tains Gibbs energies of formation for a wide variety of 
amino acids and aldehydes, in addition to water, ammo¬ 
nia, and hydrogen cyanide, all requisite species involved 
in the Strecker-type synthesis of amino acids in mete- 
oritic parent bodies. 

For each of the chemical species we are interested in 
(reactants and products), CHNOSZ contains a set of 
AG/ which change with respect to both temperature and 
pressure, spanning 0 to 1000°C and 1 to 5000 bar. We 
plot these values and perform a least squares regression. 
We model the resulting line of best fit as in Eq. 

Note that all reactant and product Gibbs energies in 
a single Strecker reaction simulation change abruptly af¬ 
ter the boiling point of water for a given pressure. Fig¬ 
ure [^illustrates this for the reactants and products of the 
Stredcer synthesis of glycine at 100 bar. In this plot, hy¬ 
drogen cyanide, formaldehyde, ammonia and glycine are 
discontinuous at 311.03°C, and also the curve represent¬ 
ing water—albeit difficult to see at this scale—slightly 
decreases in slope after this temperature. These curves 
represent a first-order phase transition; as their first 
derivatives are discontinuous. Glycine jumps to a higher 
Gibbs energy of formation after the liquid-to-gas phase 
transition of the mixture, corresponding to a decrease 
in energetic favourability. Conversely hydrogen cyanide, 
formaldehyde and ammonia become more favourable af¬ 
ter the liquid-to-gas transition due to their discontinuous 
decrease in Gibbs energies. 

Figure [^ shows the effect of varying pressure on the 
Gibbs energy of glycine. Note how each point of dis¬ 
continuity is always at the boiling point of water for the 


respective curve’s pressure. There is also an obvious lack 
of dependence on pressure before each curve reaches its 
liquid-to-gas phase transition; three pressure curves are 
shown (1.01325 bar, 50 bar, 100 bar), and they vary by 
approximately half a kilojoule from 0-100°C. (Pressure 
hardly affects aqueous solution Gibbs free energies due 
to the fact that water is nearly an incompressible fluid.) 
Therefore, in our computations, we set the pressure to 
a static value of 100 bar in order to keep the theoretical 
mixture in the aqueous phase for a larger temperature 
range (as the aqueous phase is a requirement of Strecker 
synthesis). This pressure value is on the higher end of 
interiors of planetesimals, and is in good a greement with 


Schulte & Shock 

(1995 

I, 

Cohen & Coker 

{mo 

I, and 

Warren ( 

WP2 

). in Appendix 1, we show the Gib 

□s free 

energy depeno 

encies on pressure and temperature for the 


other nine simulated amino acids (and a-aminobutyric 
acid). 



Figure 3. An example plot showing Gibbs energy dependence on 
pressure and temperature for the amino acid glycine. Temperature 
is allowed to vary from 0 °G to 500 ° C. Three curves are shown, 
representing the Gibbs energies at pressures of 1.01325 bar (blue 
curve), 50 bar (green curve), and 100 bar (red curve). Note the 
discontinuities at the boiling point of water for each pressure. 


Our thermodynamic modeling software (see Section 
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4.2 1 , requires the set of coefficients a through f for each 
chemical species, as from Eq. CHNOSZ contains the 
data necessary to calculate the coefficients for most of 
the species in our system, with the exception of glyco- 
laldehyde. We have data for all other species used in our 
model: initial reactants (water, hydrogen cyanide, and 
ammonia), five aldehydes (formaldehyde, acetaldehyde, 
butyraldehyde, propionaldehyde, and pentanaldehyde), 
and the 17 amino acids listed in Table Glycolaldehyde 
was model ed as a mixture of aceta l dehyd e and acetic 
acid, as in Espinosa-Garcfa & Dobe (2005). We obtain 
the associated G ibbs coeffic i ents to r glycolaldehyde using 


the technique in Emberson (2010). 


It is important to note t 


rat the GHNOSZ database 


does not contain Gibbs formation energies for all the 
aldehydes listed in Table [l] GHNOSZ only contains 
Gibbs data for 10 aldehydes, all differing by one carbon 
atom (the first, formaldehyde with one carbon atom, and 
the last, decanal with 10 carbon atoms). We have specific 
coefficients for formaldehyde, acetaldehyde, and glyco¬ 
laldehyde, which are the Strecker precursors for glycine, 
alanine, and serine, respectively. 

All additional aldehydes in Table are modeled using 
the closest related aldehyde for which CHNOSZ contains 
available data. The three substitute aldehydes we use are 
butyraldehyde, propionaldehyde, and pentanaldehyde. 
These three, combined with formaldehyde, acetaldehyde, 
and glycolaldehyde, provide us with a range of alde¬ 
hyde species which contain successive numbers of carbon 
atoms. The carbon count is the primary feature used to 
match Strecker aldehydes with available species substi¬ 
tutes. For example, the Strecker aldehyde necessary to 
synthesize valine is isobutyraldehyde, for which CHNOSZ 
does not contain data. We model the synthesis of valine 
using butyraldehyde rather than isobutyraldehyde, as it 
is a close substitute and contains an identical number of 
each atom. 

Aldehyde substitutions are used for val, ile, leu, lys, 
asp, glu, and thr, and are selected from the available 
aldehyde choices based on which contains the most sim¬ 
ilar molecular makeup. We require that the substitute 
aldehyde contain, at a minimum, an identical number 
of carbon atoms. We have made some quantitative es¬ 
timates on how these aldehyde substitutions affect the 
results of our Strecker synthesis reactions in Appendix 
H. Synthesis simulations for gly, ala and ser, for which 
the specific requisite aldehyde Gibbs data is used, could 
be considered the most accurate. We use butyraldehyde 
for val and pentanaldehyde for both leu and ile. The 
molecular makeups of the substitutes match identically 
the actual Strecker aldehydes in these cases (C 4 H 8 O for 
butyraldehyde, G 5 H 10 O for pentanaldehyde). The sub¬ 
stitutes aldehydes for lys, asp, glu, and thr contain the 
same number of carbon atoms as their actual aldehy¬ 
des. We use pentanaldehyde for lys, propionaldehyde 
(GaHgO) for asp and thr, and butyraldehyde for glu. Our 
analysis in Appendix H shows that our results are reason¬ 
ably accurate for all but two amino acids, glutamic acid 
and threonine, for which sensitivity to initial conditions 
could be important. 

4.2. Minimizing Gibbs Free Energy for the System: 

Thermodynamic Modeling 


We use a modeling software called ChemApp 
(distributed by GTT Technologies, http://www.gtt- 
technologies.de/newsletter) to perform a global mini¬ 
mization of Gibbs energy and calculate a theoretical yield 
of amino acids. GhemApp is a computational subroutine 
database of chemical thermodynamics which functions 
assuming a thermodynamic equilibrium of the system in 
question. For these systems, Gibbs free energy is mini¬ 
mized for a closed chemical system in equilibrium. The 
model simulates reactions in meteoritic parent bodies at 
temperatures and pressures that should theoretically ex¬ 
ist in parent body interiors as a consequence of radionu¬ 
clide decay. Thermal evolution models using radionu- 
clide decay as a heat source have also been considered by 


Grimm & McSween (1989). 

'i'he functional structure of GhemApp is to read-in a 
data file containing the relevant molecular information 
for a given trial, including thermochemical data for all 
species (coefficients a-f calculated from GHNOSZ) and 
the temperature ranges of validity. The basic compo¬ 
nents of every file are the series of basic elements in¬ 
volved in our research: carbon, hydrogen, nitrogen, and 
oxygen. Each chemical species we consider is defined as 
the addition of its elemental components. We build a 
program which reads in the data hie, and dehne a set 
of GhemApp subroutines at specihc internal conditions 
including temperature, pressure, and chemical compo¬ 
sition. In the case that reactants have far larger con¬ 
centrations than the products, it is safe to assume weak 
coupling between the reactants and do single reaction 
calculations. For this reason our ChemApp system con¬ 
sists of several individual Strecker reactions. Each indi¬ 
vidual Strecker reaction is simulated by minimizing the 
total Gibbs free energy of the system using the thermo¬ 
chemical input from each individual reaction data hie. 
Each ChemApp simulation includes all species involved 
in the specihc Strecker reaction, and is created using Fig¬ 
ure and Table As already indicated, some simula¬ 
tions include aldehyde substitutions. These substitutions 
are necessary due to the unavailability of thermochem¬ 
ical data for some aldehyde species. ChemApp’s Gibbs 
energy minimization is constrained by mass balance, and 
returns a set of molecular concentrations from the entire 
pool of reactants and products (one concentration value 
for each molecular species) which minimizes Gibbs free 
energy of the system. 

It is important to note that GhemApp has no inherent 
method for differentiating between molecular structures. 
Because each input species is broken down to its elemen¬ 
tal form, two GhemApp trials run with different molec¬ 
ular abundances, but identical elemental abundances, 
would return equivalent output concentrations. Addi¬ 
tionally, our current treatment assumes an ideal solution 
for each reaction, and we only allow the chemistry to 
be in a single phase—aqueous or gaseous—at any given 
temperature and pressure. 

Our current system restricts the possible outcomes of 
GhemApp to include only those species in which we are 
interested. Sources of initial reactants, e.g., hydrogen 
cyanide and ammonia, as well as alternate Strecker prod¬ 
ucts, are not considered. This means, for example, we do 
not include the formation pathways which lead to the for¬ 
mation of aldehydes in ChemApp; we assume some initial 
aldehyde concentration. Also, although we do not explic- 
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itly consider the hydroxy acids which Strecker synthesis 
may also produce in these calculations, w e do address 
their effects later in the paper (see Section 7.4). In this 
way, we have a limited set of initial species and a limited 
set of possible products. ChemApp output reports the 
resulting abundances of amino acids and Strecker reac¬ 
tants which minimize Gibbs energy, without considering 
alternate products which may detract from amino acid 
abundances. 

For amino acid synthesis calculations, ChemApp re¬ 
ceives as input the data file containing the Gibbs coeffi- 


i.m 

3 Gl 


for all species considered 
ibbs free energy is inde- 


cients a-f we calculated in Eq. 
in a given run. Because the 
pendent of pressure until after the liquid-to-gas phase 
transition, and we are only interested in amino acid syn¬ 
thesis via aqueous phase reactions, we set the pressure 
to a value on the upper boundary of typical parent body 
interiors—keeping the simulation in the aqueous phase 
for a wide temperature range—and hold it constant. Fac¬ 
tors which may then be varied between ChemApp runs 
are temperature and the initial concentrations of reac¬ 
tants. 


4 . 3 . Internal Conditions of Planetesimals 

Our current ChemApp treatment includes a set pres¬ 
sure at 100 bar and a temp erature range from 0°C to 
500°C in increments of 25°C. iMcSween Jr. et al. (12002 1 


and Travis & Schubert (2005) model thermal evolution 


and internal heating processes in meteoritic parent bod¬ 
ies of carbonaceous chon drites. The temperatures we 
use are in agreement with Travis & Schubert (2005) for 
the parent bodies of less altered carbonaceous chondrites. 
The higher temperatures follow the temperature c lassifi¬ 
cati o ns as l aid out in numer o us stu dies, i ncludin g Seph- 


ton (2002), Weisberg et al. (2006), and Taylor (2011), 
who discuss the temperature ranges associated with dif- 
ferent subclasses of carbonaceous chondrite. Our tem¬ 
perature ranges are also in agreement with peak alter¬ 
ation temperat ures among carbonaceou s chond rite petro¬ 
logic groups in Weiss & Elkins-Tanton (2013). 

The internal conditions tor our hducial niodel are in 
agreement with cases 1-3 of the mo deled carbonaceous 
chond rite parent bodies proposed by [Travis &: Schubert] 
(2005); that is, a spherical rocky body with a porosity of 
20%, a radius of 50 km and a rock density of 3000kg/m^; 
initially with ice water of density 917kg/m^ completely 
filling the pores of the body. 

We are primarily interested in the effects of varying 
temperature, organic composition and water on amino 
acid synthesis. As discussed in Paper I, different sub¬ 
classes of carbonaceous chondrite may be linked to dif¬ 
ferent parent body formation conditions and therefore 
different amino acid abundance levels. We vary temper¬ 
ature in our GhemApp runs to recreate the temperature 
range across meteorite classifications. 

Table lists the values and ranges of initial concentra¬ 
tions of organics input into ChemApp. These values and 
ranges were selected as a starting place based on con¬ 
centrations of organics found in comets. The molecules 
for which we have concentration data include ammo¬ 
nia, hydrogen cyanide, formaldehyde, acetaldehyde, and 
glycolaldehyde; all normalized to percent water. These 
concentrations are in agreement with cometary data in 


Table 2 

Initial concentrations of organics used in ChemApp 
model. Concentrations reported in p ercent normalized to 
water. Values are in agreement with IBocke l ee-Mo rvan et I 

ah] ll 2000 ll. lEhrenfreund fc Charn!evri[ 2 nDn _ 

[hjockelee-Mor van et ah] ||2UU4 | ). j urov i sier et ah] 20041, 
and |Mumma! & Cnarnley] l |2Ull| l. 


Molecular Species Concentration ([X]/[H 20 ]) 


Water 

Ammonia 

Hydrogen cyanide 

Formaldehyde 

Acetaldehyde 

Glycolaldehyde 

Butyr aldehyde 

Propionaldehyde 

Pentanaldehyde 

Estimated ranges. 


100 

0.7 

0.25 

0.066 

0.02-0.08 

0.0005-0.04 

0.0005-0.04* 

0.0005-0.04* 

0.0005-0.04* 


Bockelee-Morvan et al.|([2000[), Ehrenfreund & Gharnlev 

(5DD1T 

), Bockelee-Morvan et al. (I 


Grovisier et al. 

(2004 

), and 

Mumma & Gharnley 

■pm 

), which include 


Gometary concentrations were chosen as a starting 
point for the reactants within our model planetesimal 
instead of carbonaceous chondrite concentrations for 
several reasons. Comets are thought to be the mos t 


unmodified bodies in the solar system (Rauer 2008) 
and molecules in comets could have also been available 
to planetesimals at the time of the la tter’s formation 
( Schulte fc Shock||2004 Alexand^|2011 ). Carbonaceous 
chondrite matrices on the other hand are not thought to 
be pristine, being d epleted in volatiles to varying degrees 
(Bland et al. 2005). Carbonaceous chondrites are also 
from parents bodie s that have undergone h igh degrees of 
aqueous alteration ( Cobb fc Pudritz|2014 ), and are for a 
number of reasons mo re susceptible to weathering than 
other meteorite types ([Bland et al.|2006|) ; making it diffi¬ 
cult to obtain accurate parent body molecule fractions to 
water. Therefore it is safe to assume that cometary con¬ 
centrations may be more likely to represent the molecular 
concentrations in a planetesimal during the formation of 
the solar system. 

There are a few additional aldehydes included in our 
model for which we have no available cometary data. 
These are butyraldehyde, propionaldehyde, and pen¬ 
tanaldehyde. In these cases, we match their initial con¬ 
centrations with the ranges of initial concentrations for 
glycolaldehyde. We deem this to be an acceptable esti¬ 
mate because of the greater complexities of these alde¬ 
hydes compared to glycolaldehyde. 

Our first model treats a theoretical parent body as 
chemically homogeneous; the initial concentrations of re¬ 
actants are uniform throughout the model parent body. 


5 . RESULTS - TOTAL ABUNDANCES 
5 . 1 . Dependence on Temperature 

We use the physical model described above to compute 
the total theoretical yield of amino acids at temperature 
and pressure conditions analogous to those in the inte¬ 
rior of a meteoritic parent body. The results of different 
models are shown in Figured 
In Figure data points from each model are shown 
by markers located at intervals of 25°G for the tempera- 
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Figure 4. Theoretical total yield of amino acids in ppb, for a variety of models. Each data point on a line represents a modeled total 
amino acid abundance at a given temperature, in intervals of 25 °C. Amino acids included in each total concentration value are gly, ala, 
ser, asp, glu, val, thr, leu, ile, and lys. The shaded boxes indicate the observed range in amino acid abundance for each carbonaceous 
chondrite subcl ass (CIl, CMl/CRl, CM 2 . CR 2 . CV 3 /C 03 ). Dark blue s haded regions indica t e sub class classification temperatures in 
agreement withIZolenskv et al.l || 1993 ||: I Huss et al.lll 2006 ll;IKrot et a l.l(l 2006 t:IGrimm fc McSweenI lll 989 t. The aqua shaded region indicates 
cla ssification temp eratures suggested in|Morlok &: hi 5 ourel| (| 2 Ui 3 |l. Ked shaded regions correspond to temperature boundaries suggested 
by|Sephton| 1200 ^. 


ture range 0°C to 500°C and joined by a connecting line. 
Each data point represents the theoretical total amino 
acid yield. At each temperature point for each model, 
ChemApp returns some concentration value for individ¬ 
ual amino acids; we plot the total amino acid content. 
The modeled total amino acid yield at each temperature 
is a sum of the following amino acid abundances: gly, 
ala, ser, asp, glu, val, thr, leu, ile, and lys. 

The different model curves include one for a parent 
body of CM-like composition, and one for a parent body 
of CR-like composition. In two other models, also with 
CM-like composition, we scale down the water—and 
therefore the corresponding organics—by one half and 
one twentieth respectively. Then there are another three 
models with CR-like composition which correspond to 
Ex. 1, Ex. 2, and Ex. 3, in which we alter the rela¬ 
tive abundances of initial aldehydes from the CR column 
of Table Ex. 1 was run with aldehyde concentra¬ 
tions increased by an order of magnitude. Ex. 2 with 
aldehyde concentrations decreased by an order of mag¬ 
nitude, and Ex. 3 with a disrupted aldehyde pattern: 
formaldehyde, acetaldehyde, and butyraldehyde concen¬ 
trations are equal to the values shown in the CR col¬ 
umn of Table and glycolaldehyde, propionaldehyde, 
and pentanaldenyde concentrations have been increased 
by an order of magnitude. 

The observed amino acid abundances are shown as 
shaded regions in Eigure|^ The upper and lower bound¬ 
aries of each shaded region enclose the range of observed 
abundances for a carbonaceous chondrite subclass (CIl, 
CMl/CRl, CR2, CM2, CV3/C03). The left and right 
boundaries of each shaded region are the suggested tem¬ 
perature classifications for the corresponding subclass. 
Blue shaded regions cover temperature classifications in 


good a CTeement with the w ork done by Zolenskyet_^h 
( 1993 1; jjj/£g6tal (2006); Krot et al. ( |2006| ); Grimm 


& iVlcSweenllOHO). 'The aqua shaded region indicates 


the broa.d class ihcation boundaries cited by |Morlok ^ 


Libourel (2013), in which aqueous alteration may occur 
in hydrated regions between temperatures of 25° C and 
300° C. Red shaded regi ons show subcla ss classification 
boundaries suggested by Sephton (2002). 


100.00"C 

263.97"C 


— 1.01325 bar 

- - 50 bar 

100 bar 

311.03°C 


100 200 300 400 500 

Temparature (°C) 


Figure 5. Theoretical total abundances of amino acids for a par¬ 
ent body of CR-like composition, simulated at different pressures. 
The blue, green and red curves represent the total abundances 
at 1.01325 bar, 50 bar and 100 bar respectively. The three la¬ 
beled temperatures ( 100 . 00 °C, 263 . 97 °C and 311 . 03 °C) represent 
the boiling point of water at each of these three pressures. 


We observe good agreement between our CR2-type 
thermodynamic model and observed amino acid abun¬ 
dances for CR2 meteorites. In fact, all models except Ex. 
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Table 3 

Concentrations of aldehydes used in CM-type and CR-type relative frequency models which minimize and three 
experimental models of aldehyde concentrations (Ex. 1, 2 and 3). 


Molecular Species 

CM ([X]/[H20]) 

CR ([X]/[H2 0]) 

Ex. 1 ([X]/[H20]) 

Ex. 2 ([X]/[H 2 0]) 

Ex. 3 ([X]/[H 2 0]) 

formaldehyde 

0.0660 

0.0660 

0.6600 

0.0066 

0.0660 

acetaldehyde 

0.0350 

0.0726 

0.7260 

0.0073 

0.0726 

butyraldehyde 

0.0073 

0.0114 

0.1140 

0.0011 

0.0114 

glycolaldehyde 

0.0112 

0.0066 

0.0660 

0.0007 

0.0660 

propionaldehyde 

0.0042 

0.0011 

0.0110 

0.0001 

0.0110 

pentanaldehyde 

0.0066 

0.0021 

0.0210 

0.0002 

0.0210 


1 lay within the range of observed amino acid abundances 
for CR2 meteorites. We observe an approximate order 
of magnitude difference between our CM2-type thermo¬ 
dynamic model and observed amino acid abundances for 
CM2 meteorites. The total amino acid abundances for 
the CR2- and CM2-type models are practically equiva¬ 
lent. This reveals a lack of sensitivity in the resultant 
total amino acid abundances to slight variations in the 
initial aldehyde concentrations. Theoretical abundances 
remain stable as temperature increases from 0°C until 
about 200° C—where the curves begin to decrease very 
slowly until they reach the drop-off point at the liquid- 
to-gas phase transition. 

Figure illustrates how the pressure of the thermo¬ 
dynamic system plays a direct role in the temperature 
ranges at which aqueous alteration can occur. Aqueous 
alteration in petrologic types 1-3 occu rs up to 300°C, as 
suggested by Morlok & Libourel (20131. Since the liquid- 
to-gas phase change turns oft ammo acid synthesis—as 
made clear by the abrupt vertical drop-off of each curve 
in Figure at their respective boiling points of water— 
the sugg e sted l imit of aqueous alteration by [Morlok ^ 
jLibourelj (|2013|) may be a product of pressures not ex¬ 
ceeding 10ir"5ar wi thin most met eorite parent bodies. 
On the other hand, Sephton (2002) suggests that petro¬ 
logic type 2 meteorites can synthesize amino acids up 
to 400° C. In order for our simulation to yield non-zero 
abundances of amino acids at this temperature, the ini¬ 
tial pressure condition necessary would need to surpass 
the estimated pressure near the core of 10 Hygiea (the 
4th largest as teroid in the S olar System) of 246 bar; as 
calculated by Bardon (2013). 


5.2. Dependence on Water Content 

A possible reason why our CM2-type thermodynamic 
model returned abundances an approximate order of 
magnitude greater than observed amino acid abundances 
in CM2 meteorites, is that different subclasses of mete¬ 
orite could originate from different parent bodies, some of 
which tend to have less water. To exemplify these scenar¬ 
ios, the 0.5 and 0.05 water CM-type models in Figured 
calculate total amino acid abundances with one half and 
one twentieth the water content in their planetesimals 
with respect to the other models. The 0.5 water CM- 
type model decreases in total amino acid abundance by 
approximately a factor of 3 compared to the CM-type 
model, but still does not reach the boundary of the ob¬ 
served CM2 meteoritic abundances. The 0.05 water CM- 
type model on the other hand decreases in total abun¬ 
dance by over an order of magnitude compared to the 
CM-type model and is well within the range of observed 


CM2 meteoritic abundances. 

ChemApp is only able to assume a uniform composi¬ 
tion. Running separate ChemApp simulations with dif¬ 
ferent levels of water and organics—pertaining to each 
meteorite subclass in our study—might help recover the 
overall observed amino acid abundance trends we see 
among the various carbonaceous chondrite subclasses. 

6. RESULTS - RELATIVE FREQUENCIES 
6 .1. CM2 Meteorites 

Figure]^ shows the change in relative frequencies of 10 
amino acids with temperature in the CM2 subclass. Rel¬ 
ative frequencies are shown in moles of amino acid/moles 
glycine. The amino acids have been arranged in order 
of decreasing observed frequency. Each panel shows (in 
brown with black horizontal stripes) observed frequencies 
of amino acids relative to glycine among CM2 meteorites. 
Beside these observed frequencies stand (each in an indi¬ 
vidual color) the simulated relative frequencies of amino 
acids with respect to glycine. Each subsequent panel 
depicts the progression of simulated relative frequencies 
with temperature. We note that 150°C was selected as 
an arbitrary initial temperature, representative of petro¬ 
logic type 2 meteorites. As discussed in the previous 
section, our model does not recreate the low abundance 
pattern at cooler temperatures, and therefore cannot dif¬ 
ferentiate between amino acid abundances from 0°C to 
150°C. 

The ranges of initial concentrations of aldehydes go¬ 
ing into the model were taken directly from—or esti¬ 
mated using—observations of cometary abundances, as 
discussed above (see Table [^. In an effort to con¬ 
strain formation conditions within a planetesimal, we 
performed a chi-squared minimization to optimize the 
concentrations of aldehydes down to single values. The 
concentrations producing the lowest chi-squared value 
(to four decimal places), used in the model producing 
Figure [6 are shown in the CM column of Table 

The Miects of this optimization are apparent in the ex¬ 
cellent agreement between observed frequencies and the 
model frequencies at 150°C. All frequencies also remain 
well within the observed frequency error bars at 200 °C, 
but as temperature increases further, we see a more dra¬ 
matic change the frequency pattern. Threonine (thr) is 
the first amino acid to break the pattern, followed quickly 
by lysine (lys). Both amino acids produce concentrations 
less than 3*10“^ that of glycine by 300°C. At this tem¬ 
perature, isoleucine (ile) is the only other amino acid 
to have a noticeable frequency degradation. All other 
amino acids retain modeled abundances akin to their ob¬ 
served frequencies throughout the temperature progres- 
























10 


Cobb et al. 




10 ^ 


300 c 


10 '" 


10 '^ 


Gly Ala Glu Ser Asp Lys* Val Thr He Leu 



Observed 

Simulated 



Figure 6. Relative frequencies of amino acids in the CM2 subclass. Abundances of 10 amino acids are shown relative to glycine, 
reported in moles of amino acid/moles of glycine. Observed frequencies are shown in brown with black horizontal stripes, and black error 
bars representing the standard deviation. The modeled frequencies are placed beside the observed frequencies with a unique solid color 
representing each amino acid. The four panels represent the progression of relative amino acid abundances with temperature, starting at 
150°C and increasing by 50°C until 300°C. *Limited meteoric data was available for averaging the observed frequency. 


sion. This can be partially explained by thermodynamic 
favourability, as lys and ile have the highest Gibbs en¬ 
ergies of formation of the 10 amino acids through tem¬ 
perature increase, and are therefore the least favourable 
(see Appendix I Figure [l0|). 

6 .2. CR2 Meteorites 

Figure shows the change in relative frequencies of 
nine amino acids with temperature among the CR2 me¬ 
teorites. Relative frequencies are given in moles of amino 
acid/moles glycine. This time, each panel shows (in 
brown with black horizontal stripes) observed frequen¬ 
cies of amino acids relative to glycine among CR2 me¬ 
teorites. Beside these observed frequencies stand (each 
in an individual color) the simulated relative frequen¬ 
cies of amino acids with respect to glycine. The amino 
acid lysine, shown in Figure is not shown for the CR2 
equivalent, as there are no data among CR meteorites 
for lysine concentrations. 

As for the CM2 subclass, we perform a chi-squared 
minimization on aldehyde abundances for the CR2 sub¬ 
class. This is done by varying the initial aldehyde con¬ 


centrations and measuring the chi-squared between the 
relative amino acid frequencies from simulation and the 
average relative amino acid frequencies in CR2 mete¬ 
orites. The aldehyde abundances—in percent normalized 
to water—which minimize the chi-squared are shown in 
the CR column of Table |3l 

These optimized aldehyde abundances produce model 
relative frequencies in good agreement with observed rel¬ 
ative frequencies at 150°C (see Figure [^. All model fre¬ 
quencies stand well within their observe frequency error 
bars. The pattern then changes rapidly with increasing 
temperature. Some relative abundances (glu, val, leu 
and ser) increase between 150°C and 200°C, while oth¬ 
ers (ala, ile and thr) decrease. Ile and thr have the largest 
decrease in relative abundance from 150°C to 200°C. Thr 
then stabilizes whereas ile decays to 10^ by 250°C; and 
to only trace amounts by 300°C. Asp, and arguably ala, 
are the only amino acids to retain abundances similar 
to their observed frequencies throughout the shown tem¬ 
perature range. 

As these figures show frequencies relative to some base¬ 
line (in this case, glycine), it is entirely reasonable for 
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150 C 


200 C 



Gly Ala Glu Val Leu Ser Asp He* Thr* 


Gly Ala Glu Val Leu Ser Asp He* Thr* 


Figure 7. Relative frequencies of amino acids in the CR2 subclass. Abundances of nine amino acids are shown relative to glycine, 
reported in moles of amino acid/moles of glycine. Observed frequencies are shown in brown with black horizontal stripes, and black error 
bars representing the standard deviation. The modeled frequencies are placed beside the observed frequencies with a unique solid color 
representing each amino acid. The four panels represent the progression of relative amino acid abundances with temperature, starting at 
150°C and increasing by 50°C until 300°C. *Hydrolyzed amino acid abundances were not determined due to co-elutions, therefore the 
relative abundances represent those obtained from the unhydrolyzed samples (moles of free amino acid/moles of free glycine). 


future calculations to scale overall abundances of alde¬ 
hydes. We assume some base concentration of formalde¬ 
hyde, as this is the precursor for glycine, to which we 
normalize additional amino acids. It is the abundances of 
the aldehydes relative to formaldehyde that are respon¬ 
sible for the agreement between observed and modeled 
amino acid frequencies. 

In the case of Figure the abundance of acetalde¬ 
hyde is higher than its abundance in the CM2 model by 
approximately a factor of 2. This makes sense visually. 
Acetaldehyde is the aldehyde involved in the modeled 
synthesis of ala; and the CR frequency of ala is approxi¬ 
mately twice as much as its CM counterpart. 

7. DISCUSSION 
7.1. Minimization 

The aldehyde concentrations listed in Table mini¬ 
mize x^ over variations in both concentration and tem¬ 
perature. While this helps to constrain formation condi¬ 
tions of amino acids in planetesimals, we are left with no 
impression of which thermodynamic consideration has a 


larger impact on organic synthesis. There are three as¬ 
pects in our thermodynamic model that may be manip¬ 
ulated: pressure, temperature, and composition. 

We set pressure to a typical value for the interior of par¬ 
ent bodies and hold it constant. The effect of pressure 
on Gibbs free energy before the liquid-to-gas phase tran¬ 
sition is minimal in comparison to the effects of changing 
temperature, as discussed above. 

With regard to which aspect, temperature or compo¬ 
sition (which we treat, for now, as a simplified model of 
aldehyde abundances), has a more significant effect on 
relative frequencies, we consider the evolution of x^ val¬ 
ues. We consider the relative frequencies of amino acids 
in CR meteorites as a representative case. 

Figurej^shows the evolution of x^ with increasing tem¬ 
perature, for a variety of initial aldehyde concentrations. 
The separate plots represent different ChemApp models 
which were run with aldehyde concentrations listed in 
Table The plot labeled ‘CR-Type’ shows the x^ val¬ 
ues for the model we used in Figure Examples I, 2, 
and 3 show evolution for three different combinations 
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Figure 8. values for different ChemApp models, reflecting goodness-of-fit evolution for temperatures ranging from 150°C to 300°C, 
for four different initial aldehyde concentrations. The CR-Type plot represents the x^ values for the aldehyde abundances shown in the CR 
column of Table and used in FigureThe Ex. 1 plot represents the x^ values for a model run with aldehyde concentrations an order of 
magnitude greater than those shown m the CR column of Table ^ The Ex. 2 plot represents the x^ values for a model run with aldehyde 
concentrations an order of magnitude lower than those shown in the CR column of Table The Ex. 3 plot represents the x^ values for a 
model run with the formaldehyde, acetaldehyde and butryaldehyde concentrations shown in the CR column of Table|^ and glycolaldehyde, 
propionaldehyde and pentanaldehyde concentrations greater by an order of magnitude. The red dotted line represents the x^ value (from 
the x^ distribution table) for a p-value of 0.05 with 7 degrees of freedom. Any x^ values lying above the dotted line represent a rejected 
null hypothesis. 


of aldehyde concentrations, all varying with respect to 
the CR-Type concentrations. The concentrations used 
in ‘Ex. 1’ are one order of magnitude greater than those 
in the CR column of Table and those used in ‘Ex. 2’ 
are decreased by one order ofmagnitude. ‘Ex. 3’ disrupts 
the relative abundance pattern; concentrations of three 
aldehydes (formaldehyde, acetaldehyde, butyraldehyde) 
are those in the CR column of Table but the three 
least abundant aldehydes (glycolaldeh 3 ^e, propionalde¬ 
hyde, pentanaldehyde) have been increased by one order 
of magnitude. 

The red dotted line represents a value of 14.067. 
This value comes from the distribution table for a 
p-value of 0.05, with n - 1 = 7 degrees of freedom, (n 
represents the number of relative amino acid frequencies 
which the simulation can vary; see Figure [H) All 
values higher than this line represent a rejected null hy¬ 
pothesis, i.e., a statistical unlikelihood of obtaining those 
simulated relative abundances given the observed values. 

The result of this simple x^ comparison shows that ini¬ 
tial aldehyde composition has a greater effect on amino 
acid frequencies than does temperature evolution. In the 
CR-Type and Ex. 2 cases, x^ increases by approximately 
a factor of 3 over a 200°C change. In the case of Ex. 1 
and Ex. 3, x^ decreases by approximately a factor of 2 
over a 200°C change. Compositional effects on the con¬ 


trary are usually much larger. The change in between 
the CR-Type composition and Ex. 1 and 3 compositions 
at 150°C is nearing 2 and 3 orders of magnitude respec¬ 
tively. Only the Ex. 2 composition doesn’t change in 
relative to the CR-Type composition. The similar¬ 
ity between the CR-Type and Ex. 2 curves illustrates 
that relative amino acid abundances to glycine tend not 
to change very much if you equally scale aldehyde con¬ 
centrations down (as opposed to up) from the CR-Type 
concentrations. 

Both the CR-Type and Ex. 2 curves also stay well 
below the null hypothesis rejection line—yielding statis¬ 
tically probable relative frequencies. The Ex. 1 and Ex. 
3 curves on the other hand start and remain above this 
line. This further illustrates the sensitivity of increas¬ 
ing the initial aldehyde concentrations versus increasing 
temperature. 


7.2. Does the Onion Shell Model Aecount for Amino 
Acid Data? 

As discussed in Paper I, the onion shell model for 
planetesimals proposes that the carbonaceous chondrites 
within a cert ain subclass may share the same meteorit e 
parent body (Ehrenfreund et al. 2001 Glavin et al. 2011), 
and that the different petrologic types within each sub- 
class originate from different layers within the parent 
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Simulation 
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CR2 observed data 
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Figure 9. Theoretical total yield of amino acids in ppb with varying planetesimal water content. Water content varies in weight percentage 
over the entire planetesimal from 0% to 50%. Simulation is run with CR-like initial concentrations at lOObar and 100°C. The shaded boxes 
indicate the observed range in amino acid abundance for each carbonaceous chondrite subclass (Cll, CMl/CRl, CM2, CR2, CV3/C03). 
The individual data points represent the percent uncombined wate r co ntent of sepa rate samples of some of the carbonaceous chondrites 
used in this study, measured and/or reported by | Jarosewich| (|1990|| and|Wiik| (|1956|l. The percent uncombined water content of the data 
points in the left panel are as stated in the literature. The percent uncombined water content of the data points in the right panel are 
scaled down by a factor of 100. 


body; wit h each layer corresponding to a different tem- 


20131. 


perature ( Wejss fc Elkins-Tanton 
In Figure the change in theoretical total amino 
acid abundances over the aqueous phase temperature 
range is too minute to support a simple temperature- 
differentiated onion shell model. For example the CRl 
chondrites have total abundances just above 10^ Ppb. 
If these meteorites are assumed to originate from loca¬ 
tions within a parent body relatively closer to the core 
(whereas the CR2 meteorites, with total abundances near 
10 ^ ppb, would be in the in-between layers associated 
with cooler temperatures), we would expect our CR2 
models in Figurej^to move towards 10^ ppb as tempera¬ 
tures increase. Similarly, the CR3 chondrites, associated 
with the coolest alteration temperatures (pertaining to 
locations within the parent body nearest to the surface), 
also have a slightly lower abundance than the CR2 me¬ 
teorites (near 10^ ppb for the only analyzed CR3 mete¬ 
orite, QUE 99177). However, our CR2 model shown in 
Figure H] does not tend towards lower abundances at the 
highest temperatures (near 300° C), or at the lowest tem¬ 
peratures (near 25°C) thought to be in meteorite pareii t 
bodies of petrologic types 1-3 (|Morlok & Libourel|2013|). 
The lack of such behavior suggests that one of the mam 
drivers of amino acid synthesis in planetesimals is the 
amount of water and organics in the planetesimal and 
not the temperatures inside each layer. 


7.3. Water Content and Amino Acid Abundance 

Figure explores the effect of varying water content 
further by showing the dependence of theoretical to¬ 
tal amino acid yields to overall water percentage in the 
model planetesimal. The simulation is done with CR-like 


aldehyde concentrations, at 100 bar and 100°C. These 
initial conditions are fairly arbitrary as Figure shows 
how abundances are relatively the same across aque¬ 
ous phase temperatures, and CR-like and CM-like initial 
compositions also produce similar curves. 

The black curve in both panels of Figure [^illustrates 
how total amino acid abundance can changeoy varying 
the water content of the fiducial model planetesimal. The 
shaded boxes represent the minimum and maximum ob¬ 
served abundances in each of the carbonaceous chondrite 
subclasses, fit onto the simulation line. By drawing a line 
downwards from the bottom-left and top-right corners of 
each box, one could determine the theoretical minimum 
and maximum percent water content in the parent bodies 
hosting each of the carbonaceous chondrite subclasses. 
By scaling the composition of a planetesimal from 0% 
water and 100% rock to 50% water and 50% rock, the 
simulated total amino acid abundances can be fit into the 
ranges of observed abundances for all the carbonaceous 
chondrite subclasses in this study. This result supports a 
composition-differentiated model; specifically, planetesi¬ 
mals with different overall water concentrations. 

The aqua, dark blue, and red data points correspond 
to percent uncombine d w ater m e asure d and/or reported 
by Jarosewich (1990) and jWiik (1956) in separate sam- 


ples (but the same body) of soine of the carbonaceous 
chondrites that have been used in this study. Since the 
amino acid abundances and uncombined water content 
were not measured on the same sample, the amino acid 
abundance values for each data point correspond to the 
samples measured from the sources outlined in Paper I; 
with error bars corresponding to two or more sample 
measurements. The percent uncombined water values 
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in Jarosewich (1990) and Wiik (1956) are measured by 
heating an uncovered sarnpleto lUr C for 1 hour and 
measuring the loss in weight. The measured values for 
these data points therefore consist of both excess free 
water from weathering and original free water. 

Notice in the left panel of Figure P how all of the 
data points lie to the right of the simulation line. This 
overabundance in measured uncombined water percent 
in comparison to our simul ated values may be the result 
of substantial weathering. [Bland et al. (2006) conclude 
that for a number of reasons, carbonaceous chondrites 
are more susceptible to weathering than other meteorite 
types and that even observed falls (as opposed to finds) 
may display the effects of significant weathering. The 
timescale for atmospheric, ground and/or rain water to 
be ta ken in by rece ntly fallen meteorites could be very 
short: |Pisani| (|1864D dehydrated a sample of the Orgueil 
meteorite used in our study, and found that after only 
a few hours it had regained almost all the hygroscopic 
water that it had lost (9.15 wt%). 

In attempts to quantify the amount of weathering en¬ 
dured by each meteorite—and thus the corresponding 
excess free water in each sample—we look at the ratio 
of oxidized iron (Fe 203 ) to metallic iron (Fe(m)). The 
oxidation of Fe(m) into Fe 20 3 is a well known weath¬ 
ering effect ( Jarosewidi||1990 ), therefore this ratio could 
provide a relative scaling factor for each of the water con¬ 
tent values of our data points. Unfortunately this idea 
fails as none of our eight data samples yield non-zero and 
non-trace values for both oxidized and metallic iron. 

In the right panel of Figurej^we apply a constant scal¬ 
ing factor of 0.01 to all the uncombined water percent val¬ 
ues. This allows us to estimate the ratio of original free 
water to excess free water from weathering that would 
allow the data points to coincide better with our simula¬ 
tion curve. Since individual samples will be subjected to 
distinct amounts of weathering, this 0.01 scaling factor 
is a simplification in attempts to typify the magnitude 
of excess water due to weathering in carbonaceous chon¬ 
drites. 

7.4. Chemical Competition Between Amino Acids and 
Hydroxy Acids 

So far we have, for simplicity, focused on the forma¬ 
tion of amino acids without discussing other competing 
chemical pathways. This cannot be ignored in general, 
and we now turn to a discussion of one of the main com¬ 
peting processes. Hydroxy acids are similar to amino 
acids, except in place of the amino group (NH 2 ), they 
have a hydroxyl group (OH). The hydroxy acid analogous 
to glycine is glycolic acid, and the hydroxy acid analo¬ 
gous to alanine is lactic acid Since hydroxy a cids can 
also form from the Strecker synthesis pathway (Schulte 


&: Shock 
mation o 


1995), it is important to consider how the tor- 
hydroxy acids can reduce amino acid synthesis 


within meteorite parent bodies. A key point about this 
chemistry is the role of ammonia that is so essential for 
Strecker synthesis. As NH 3 concentrations decrease, one 
expects that amino acid p roduction will fall with r espect 
to hydroxy acid synthesis |Schulte & Shock] (|1995D. 

Table lists the hydroxy acid to analogous ammo acid 
meteorite abundance fractions for these molecules, from 
five meteorite samples. The molecula r abundances are 
obtained from Cobb & Pudritz (2014), Pizzarello et al. 


(2010) and Monroe & Pizzarello (2011). There is a large 
variation ih the hydroxy acid to analogous amino acid 
fractions between the CIl, CM2 and CR2 meteorite sub¬ 
classes. Within CR2 meteorites, there is only 0.04-0.09 
times as many hydroxy acids as their analogous amino 
acids. In CM2 meteorites, hydroxy acids are 1-28 times 
more abundant than their analogous amino acids. And 
most signihcantly, in the Ivuna CIl meteorite, hydroxy 
acids are 121-274 times as abundant as their analogous 
amino acids. 

Within CIl meteorite parent bodies, the molecu¬ 
lar fractions suggest that conditions could be more 
favourable for hydroxy acid formation over amino acid 
formation. If hydroxy acid production is favoured within 
CII meteorite parent bodies, it could explain why our 
simulations produce almost 3 orders of magnitude more 
total amino acids than are present in CIl meteorites. 
This suggests that observed NH 3 concentrations are pos¬ 
itively correlated with amino acid synthesis within CIl 
meteorite parent bodies, which is to be expected for 
Strecker reactions. 

Within CR2 meteorite parent bodies on the other 
hand, there appears to be only a small fraction of hy¬ 
droxy acids relative to their analogous amino acids. This 
could mean that there is sufficient NH 3 within CR2 par¬ 
ent bodies allowing for optimal amino acid production. 
Any variation in amino acid abundances between the 
various CR2 meteorites could then be mainly from vari¬ 
ations in o vera ll parent body water content (as shown 
in Section |7.3 ). This is validated by our simulations. 


which produce total amino acid abundances that agree 
well with the CR2 meteoritic abundances. 

The CM2 meteorites may represent the intermediate 
case between these two extremes. Here, hydroxy acid to 
analogous amino acid fractions are in between the values 
for the CIl and CR2 meteorites. NH 3 concentrations 
could therefore be moderately affecting amino acid syn¬ 
thesis within CM2 meteorite parent bodies, leading to 
similar abundances of hydroxy and amino acids. This 
fits in well with the story thus far, as our simulations 
produce about an order of magnitude more amino acids 
than the CM2 meteoritic abundances, which could be ex¬ 
plained by a moderate reduction of amino acid synthesis 
from relatively low NH 3 concentrations. 

To investigate the possible role of reduced NH 3 as the 
ultimate cause of reducing amino acid synthesis within 
CIl and CM2 meteorite parent bodies, we collected the 
available data on the ammonia content for each of the 
five meteorite sa mples in Table These ab und ances 
are obtained from Monroe & Pizzarello (2011) and Piz¬ 


zarello & Holmes (2UU9r Results show that the CR 2 
meteorites have over an order of magnitude more NH 3 
than the CM2 meteorites, and about a factor of 3 more 
NH 3 than the CIl meteorite. Since CR2 meteorite par¬ 
ent bodies have the greatest amount of NH 3 , amino acid 
synthesis should be favo ured over hydroxy acid synthesis 
( Schulte fc Shock||l995 ), which aligns well with our sim- 
ulations and the meteoritic record. The CM2 meteorites 
however—especially in the case of the Bells sample—have 
less NH 3 than the CIl meteorite. This does not fit per¬ 
fectly with this picture, where amino acid abundances in 
CIl meteorite parent bodies would seem to have been re¬ 
duced more than in CM2 meteorite parent bodies, judg¬ 
ing by the hydroxy acid to analogous amino acid ratios. 
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Table 4 

Molecular fractions of hydroxy acids to their analogous amino acids and the amm onia content within one 
CIl meteo rite, two CM2 meteorit e s, and two CR2 meteorites. Dat a obtained fromICobb fc Pud rit^||2014||. 
|PizzareIlo et al.1 | |2010| l , [Monroe fc Pizzarello| | |2011| | and jPizzarello &: Holmesj pUUUj l. 


Meteorite 

Subclass 

[mol glycolic acid]/[mol glycine] 

[mol lactic acid]/[mol alanine] 

NH 3 (nmol/g) 

Ivuna 

CIl 

120.88 

274.17 

5300 

Bells 

CM2 

16.06 

28.00 

280 

Murchison 

CM2 

0.68 

1.94 

1100 

GRA 95229 

CR2 

0.09 

0.05 

18850 

LAP 02342 

CR2 

0.04 

0.04 

14080 


This could perhaps partially be explained by the fact that 
water content within the petrographic type 1 meteorites 
tends to be higher than within the petrographic type 2 
meteorites, which could have an effect on how diluted the 
ammonia is within the CIl parent body. 


7.5. Amino Acid Decomposition 

One limitation of our model is that amino acid decom¬ 
position pathways are not simulated. Amino acids deam- 
inate and decarboxylate in aqueous solutions at rates 
that vary with temper ature ( McCollom|2013 Sato et al. 
'2004 Qian et al.|1993). High-pressure, high-temperature 


laboratory experiments infer through extrapolation that 
for temperatures less than ~60°C, some amino acids have 
half lives > 10® years. These experiments also show that 
amino acids degrade rapidly at temperatures > 200° C. 
However, since at lower temperatures, amino acid decom¬ 
position reaction rates may be longer than the a queous 
lifetime of a meteo rite parent body (~10® years) (Travis 


& Schubert 2005), it is reasonable to assume that not 
too much decay has occurred. The main exception would 
be if a meteorite parent body interior reaches termera- 
tures > 200°C, but as we have shown in Figures IQ and 
simulated relative frequency patterns diverge from the 
meteoritic values at 250° C. Hence, since we expect that 
interior temperatures never reach above 200° C, decay ef¬ 
fects should be minimal. 


7.6. A New Picture - Linking Amino Acid Synthesis in 
Planetesimals to Disk Chemistry 

The results of our simulations suggest a new picture of 
amino acid synthesis within planetesimals. Total amino 
acid abundances within planetesimals remain relatively 
stable throughout the aqueous phase temperature range. 
This suggests that temperature and pressure are not the 
main driving factors for differences in amino acid abun¬ 
dances across meteorite types. We have already argued 
that hydroxy acid production could be favoured over 
amino acid production within certain meteorite parent 
bodies (e.g. CIl and CM2), while amino acid production 
could be favoured over hydroxy acid production within 
others (e.g. CR2). This could be due to varying NH 3 to 
water levels within these parent bodies. Secondly, initial 
water content within the planetesimal, and the scaling or¬ 
ganics within that water, could account for the variation 
of amino acid synthesis within meteorite parent bodies of 
the same subclass. This is demonstrated by the shaded 
boxes in Figure which show that varying the wa¬ 
ter content within CR2 meteorite parent bodies—whose 
amino acid abundances may not be curbed by low NH 3 
concentrations—from 0.02-18 wt% (purple box), could 
account for the variation in CR2 meteoritic amino acid 


abundances. Amino acid synthesis within CIl, CMl, 
CRl, COS and CVS meteorite parent bodies on the other 
hand, could first be reduced due to relatively low NH 3 
concentrations, and then varied further by the ranges in 
their parent body water contents. 

The variation of water and ammonia within planetes¬ 
imals is ultimately tied to where and how these bodies 
formed in their host protostellar disks. Astrochemistry 
models for such systems are starting to become sophisti- 
cated enough to allow us to explore this connection (see 


Pontoppidan et al. (2014)). 

(Jur proposed picture conforms well with the water 
percen tages measured i n vari ous chondrites in the liter¬ 
ature. Weisberg et al. (2006) state that petrologic type 
1 chondrites have the highest water content by weight 
(18-22%), followed by type 2 chondrites (2-16%) and 
then type 3 chondrites (0.3-3%). Figurej^shows that in¬ 
creasing the water content in the meteorite parent body 
leads to greater amino acid production. Since CIl mete¬ 
orite parent bodies appear to have had the largest water 
percentages, the reduction of amino acid production from 
low NH 3 concentrations within these parent bodies would 
have to be greater in order to produce the comparatively 
low amino acid abundances measured in CIl meteorites. 
CR2 meteorite parent bodies on the other hand, which 
appear to have had mid-range water concentrations ( 2 - 
16%), have ranges in meteoritic amino acid abundances 
that can be explained almost completely with varying 
the water content within our simulations (from 0.2-18%). 
CMl, CRl, C03 and CV3 meteorite parent bodies would 
then require some level of amino acid abundance reduc¬ 
tion (by having relatively low NH 3 concentrations), as 
well as an assortment of water contents, to account for 
their lower and varied meteoritic am ino acid abundances . 

In circumstellar disk m odels from Woitke et al. (2009) 
and Visser et al. (2012), ice water content aiorig the 
mid-plane is found to increase swiftly (but not instan¬ 
taneously) at the ice line, after which it changes only 
minimally for several AU. The water content within each 
forming planetesimal varies as a function of its distance 
from the ice line. For the Solar System, the ice line dur¬ 
ing planet formation may have been somewhere within 
the asteroid belt, as modeled by [Raymond et al.| ([2004) 
(at 2 and 2.5 AU from the Sun)l Most of Earth s me- 
teorites (including carbonaceous chondrites) are believed 
to originate from the asteroid belt. If the water ice dis¬ 
tribution is inherited from the cicumstellar disk models 
explained above, then the proposed variation in water 
content between the parent bodies of the carbonaceous 
chondrites in our study could be possible. 

The variation of ammonia acros s protostellar disks 


has also been computed (see review Pontoppidan et al. 
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(2014)). It is known observationally that NH 3 can be 
ruled out as a significant carrier of nitrogen in the warm 
inner regions of disks so that most of the nitrogen in the 
terrestrial planet forming region may be stored as N 2 . 
At larger disk radii, nitrogen will be found in the form of 
ammonia, which explains its abundance in comets and in 
the atmospheres of giant planets. Applied to our model, 
this suggests that if high ammonia levels secure domi¬ 
nant amino acid abundances, then the parent bodies for 
such meteorites would have derived from such outer disk 
regions beyond the terrestrial planet formation region. 


8 . CONCLUSIONS 

In general, we see good agreement between our thermo¬ 
dynamic computations and observed amino acid abun¬ 
dances. When applying our fiducial model, thermody¬ 
namics well explains the total amino acid concentrations 
we see in CR2 meteorites and the relative frequencies we 
see in both the CM2 and CR2 subclasses. 

Our models for CM-type and CR-type total amino acid 
abundances at temperatures less than 311.03° C remain 
near 7x10^ parts-per-billion (ppb). This value is in agree¬ 
ment with the average observed abundance in CR2 mete¬ 
orites of 4 ± 7x10®, but approximately an order of magni¬ 
tude higher than the average observed abundance in CM2 
meteorites of 2 ± 2x10"^. These simulation results are not 
sensitive to slight changes in the initial aldehyde concen¬ 
trations. For this reason we think the results will remain 
fairly similar regardless of whether relative amino acid 
abundances in hydrolyzed or unhydrolyzed meteorite ex¬ 
tracts are used in calculations. An increase in the pro¬ 
duction of the competing hydroxy acids in the meteorite 
parent body can explain the discrepancy between the 
CM-type model computations and the observed amino 
acid abundances in CM2 meteorites. 

There is a sharp drop in total amino acid abundance 
at the liquid-to-gas phase transition of the aqueous solu¬ 
tion. We attribute this drop-off to the decline in thermo¬ 
dynamic favourability of synthesizing amino acids once 
in a gaseous phase mixture. 

Our models do not decrease in total amino acid 
abundance by any significant amount over the course 
of the aqueous phase temperature range. This lack 
of declination does not support a simple temperature- 
differentiated onion shell model. 

Instead, we propose a new explanation for amino acid 
abundances in planetesimals. In this new model, the 
production of hydroxy acids is favoured over the pro¬ 
duction of amino acids within certain meteorite parent 
bodies (e.g. CIl and CM2). This is likely due to the low 
NH 3 concentrations within these parent bodies, which, as 
we have seen, results in lower amino acid production in 
preference to hydroxy acid production. (Protoplanetary 
disk astrochemistry models predict low ammonia levels 
in the inner regions of disks.) Then, to get the variation 
in the amino acid abundances within the same carbona¬ 
ceous chondrite subclasses, the total water content (and 
scaling reactant content) could vary from parent body to 
parent body. 

Our model accounts for the relative frequency pattern 
in amino acids among CM2 and CR2 chondrites. The ob¬ 
served and theoretical frequencies with respect to glycine 
match to well within an order of magnitude in both cases, 
at a temperature of 150°C. The degradation of lysine 


and threonine for CM-type relative frequencies and of 
isoleucine for CR-type relative frequencies by 250°C and 
300° C, is in agreement with 3D models of interiors of 
young planetesimals; which are not heated much beyond 
200° C. This change in amino acid frequencies may sug¬ 
gest a method to estimate internal meteorite parent body 
temperatures within that range. 

These various aspects of the amino acid record in 
meteorites point towards the deep connection that likely 
exists with local disk astrochemistry. In particular, 
the processes that determine the abundances of am¬ 
monia, water, and aldehyde concentrations at different 
disk radii and their incorporation into planetesimals 
ultimately provide the initial conditions for Strecker 
synthesis of the amino acids. It is therefore possible 
to link biomolecule patterns in various types of me¬ 
teorites to a much broader picture of disk astrochemistry. 
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APPENDIX I - ADDITIONAL GIBBS ENERGIES OF 
FORMATION 

In Eigure we extend Figure (Gibbs free energy 
of formation for Glycine) to the other nine amino acids 
considered in our simulation, and order them by increas¬ 
ing Gibbs free energy (at fiducial values of 50°G and 100 
bar). All amino acid curves are piecewise, and increase in 
Gibbs free energy at the point of discontinuity. Table 
displays exact values for the Gibbs free energies of for¬ 
mation for each of the proteinogenic amino acids used in 
our simulation, plus one non-proteinogenic amino acid, 
at the fiducial temperature and pressure. 

It is worth noting that the order of amino acids in Ta¬ 
ble does not directly match the order of amino acids 
in Figure]^ or Figure This implies that the amount 
of each amino acid synthesized in planetesimals is dic¬ 
tated by more than just its energetic favourability. Addi¬ 
tional factors, such as the concentration of each aldehyde 
precursor, cause a dispersion in the ordering of relative 
amino acid abundances. 

q;-Aminobutryric Acid 

Figure pj] shows the Gibbs free energy of formation 
for a-aminobutryric acid with respect to temperature, 
at three different pressures. In contrast to the plots in 
Figure the Gibbs free energy of a-aminobutryric acid 
decreases after the boiling points of water for each of 
the three pressure curves. This signifies an increase in 
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Table 5 

Gibbs free energies of formation for each of the 
proteinogenic amino acids used in our simulation plus 
one non-proteinogenic amino acid; at 50° C and 100 bar. 
Amino acids are ranked by increasing Gibbs free energy 
of formation. 


Amino Acid 

Gibbs Energy (kj/mol) 

Glutamic Acid 

-730.0406 

Aspartic Acid 

-727.2654 

Serine 

-523.2076 

Threonine 

-505.6537 

Glycine 

-384.4235 

Alanine 

-375.3493 

a-Aminobutyric Acid 

-369.0258 

Valine 

-361.3631 

Leucine 

-357.0445 

Isoleucine 

-347.7252 

Lysine 

-342.9651 


Non-proteinogenic. 


energetic favourability after the liquid-to-gas phase tran¬ 
sition. 

Since we have been basing our simulations on an aque¬ 
ous phase chemical reaction (Strecker synthesis), com¬ 
paring the Gibbs energy of a-aminobutryric acid in the 
liquid phase to that of the proteinogenic amino acids 
is most important. This is shown in Table where 
a-aminobutryric acid ranks 7th in energetic favourabil¬ 
ity at 50°C and 100 bar. From compiling amino acid 
abundances from the sources outlined in Paper I, a- 
aminobutryric acid is found to be tied for 5th in highest 
average amino acid abundance in both the CM2 and CR2 
subclasses. Observed a-aminobutryric acid abundances 
therefore rank near what we would expect based on their 
relative Gibbs free energies. But as stated above, there 
are more factors to amino acid synthesis than just ener¬ 
getic favourability. 

APPENDIX II - ALDEHYDE SUBSTITUTIONS 

7 of the 10 amino acid synthesis simulations are per¬ 
formed using the Gibbs free energy data from aldehyde 
substitutes rather than the actual Strecker precursors 
listed in Table Due to the difference in steric configu¬ 
rations between the substitutes and the actual reactant 
aldehydes, we can assume that their Gibbs free energies 
may also differ. In order to try to identify the difference 
between their Gibbs free energies, we consulted the Ayers 
group in the Dept, of Chemistry and Chemical Biology 
at McMaster University, who performed an analysis us¬ 
ing an electronic structure program called Gassian (F. 
Zadeh, P. Ayers and A. Patel 2015, personal communi¬ 
cation) . 

From their analysis, it is discovered that the substi¬ 
tutes which have exactly the same atoms as the actual 
aldehydes, i.e., isobutryaldehyde (val), 3-methyl butanal 
(leu) and 2-methyl butanal (ile), have very similar Gibbs 
free energies to the actual aldehydes. Therefore the accu¬ 
racy of the val, leu and ile synthesis simulations are likely 
not greatly affected by the aldehyde substitutions. The 
substitutes which do not have exactly the same atoms 
as the actual aldehydes, i.e., 3-oxopropanoic acid (asp), 
5-aminopentanal (lys), 4-oxobutanoic acid (glu) and lac- 
taldehyde (thr), were found to have a considerable dif¬ 
ference in Gibbs free energy to the actual Strecker pre¬ 
cursors. 


In attempts to identify how greatly these four aldehyde 
substitutions might affect our amino acid synthesis re¬ 
sults, we use the rough Gibbs data computed from Zadeh 
et al. (2015, personal communication) to adjust our alde¬ 
hyde substitute Gibbs coefficients. We then rerun the 
CM-type simulations to compare to our previous results. 
We find that the reactions involving 3-oxopropanoic acid 
(asp) and 5-aminopentanal (lys) are not sensitive to the 
adjusted aldehyde Gibbs energies. 

Thus, there are only two reactions in this set, involving 
4-oxobutanoic acid (glu) and 2-hydroxy-propanal (thr), 
that are quite sensitive to adjustments in Gibbs energies. 
We therefore conclude that the aldehyde substitutions 
used for the glutamic acid and threonine simulations may 
lead to significant deviations from accuracy, whereas the 
other five aldehyde substitutions are likely acceptable. 
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Figure 10. Plots extending Figurel^for the nine other amino acids considered in this simulation, in order of increasing Gibbs free energy 
of formation (at 50°C and 100 bar). This plot depicts each amino acid’s change in Gibbs free energy of formation (kj/mol) with respect to 
temperature (°C), at various pressures. The vertical grey dotted lines represent the temperatures of the the aqueous-to-gas phase transitions 
of the amino acid (and hence of water) at 1.01325 bar, 50 bar and 100 bar—which are 100°C, 263.97°C and 311.03°C respectively. 



Figure 11. A plot of Gibbs energy dependence on pressure and temperature for the non-proteinogenic amino acid Q;-aminobutyric acid. 
Temperature is allowed to vary from 0°C to 500°C. Three curves are shown, representing the Gibbs energies at pressures of 1.01325 bar 
(blue curve), 50 bar (green curve), and 100 bar (red curve). 














































































